THE CARMICHAEL NUMBERS UP TO lO^^ 



- - - RICHARD G.E. PINCH 

O 

C !) ■ Abstract. We extend our previous computations to show that there are 1401644 

^Nl ' Carmichael numbers up to 10^^. As before, the numbers were generated 

J, ' by a back-tracking search for possible prime factorisations together with a 

O , , "large prime variation". We present further statistics on the distribution of 

Carmichael numbers. 
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1. Introduction 

\r^ , A Carmichael number N is a, composite number N with the property that for 

every b prime to N we have b^~^ = 1 mod N. It foUows that a Carmichael number 
(-H ] N must be square-free, with at least three prime factors, and that p — l\N — 1 for 

every prime p dividing N: conversely, any such N must be a Carmichael number. 

For background on Carmichael numbers and details of previous computations 
we refer to our previous paper 0]: in that paper we described the computation of 
the Carmichael numbers up to 10"'^^ and presented some statistics. These computa- 
tions have since been extended to 10""^^ |S], 10"'^'' |S] and now to 10^^, using similar 
K*" I techniques, and we present further statistics. 

1^ ■ The complete list of Carmichael numbers up to 10^^ is available from the author's 

L^ [ web site at www.chalcedon.demon.co.uk. 
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2. Organisation of the search 



P^ I We used improved versions of strategies first described in j^ . 

^^ . The principal search was a depth-first back-tracking search over possible se- 

^-^ I quences of primes factors pi,...,pd- Put Pr = IlLiPii Qr = I\i=r+iPi ^^^^ 

C^ ' Lr = 1cm {pi — 1 : « = 1, . . . ,r}. We find that Qr must satisfy the congruence 

S I N = PrQr = 1 mod Lr and so in particular Qd — Pd must satisfy a congruence 

modulo Ld-i- further pd — 1 must be a factor of Pd-i — 1- We modified this to 
terminate the search early at some level r if the modulus Lr is large enough to limit 
the possible values of Qr: which may then be factorised directly. 
H \ We also employed the variant based on proposition 2 of 0| which determines the 

finitely many possible pairs {pd-i,Pd) from Pd-2- In practice this was useful only 
when d = 3 allowing us to determine the complete list of Carmichael numbers with 
three prime factors up to 10"'^^. The count of these is 35586 (not 35585 as incorrectly 
reported in our preliminary work) . 

Finally we employed a different search over large values of pd, in the range 10^ < 
Pd < 10^, using the property that Pd-i = 1 mod (pd— 1) and factorising the numbers 
in this arithmetic progression less than IQ^^ /pd- 

3. Statistics 

We have shown that there are 1401644 Carmichael numbers up to 10^*, all with 
at most 11 prime factors. We let C{X) denote the number of Carmichael numbers 
less than X and C{d,X) denote the number with exactly d prime factors. Tabled 
gives the values of C{X) and Table El the values of C{d,X) for X in powers of 10 
up to 10^^. 
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n 


C (10") 


3 


1 


4 


7 


5 


16 


6 


43 


7 


105 


8 


255 


9 


646 


10 


1547 


11 


3605 


12 


8241 


13 


19279 


14 


44706 


15 


105212 


16 


246683 


17 


585355 


18 


1401644 



Table 1. Distribution of Carmichael numbers up to 10 
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6 
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105 
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84 


144 


27 




















255 
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172 


314 


146 


14 

















646 


10 


335 


619 


492 


99 
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1547 


11 


590 


1179 


1336 


459 


41 














3605 


12 


1000 


2102 


3156 


1714 


262 


7 











8241 


13 


1858 


3639 


7082 


5270 


1340 


89 


1 








19279 


14 


3284 


6042 


14938 


14401 


5359 


655 


27 








44706 


15 


6083 


9938 


29282 


36907 


19210 


3622 


170 








105212 


16 


10816 


16202 


55012 


86696 


60150 


16348 


1436 


23 





246683 


17 


19539 


25758 


100707 


194306 


172234 


63635 


8835 


340 


1 


585355 


18 


35586 


40685 


178063 


414660 


460553 


223997 


44993 


3058 


49 


1401644 



Table 2. Values of C{X) and C{d, X) for d < 10 and X in powers 
of 10 up to 10^^ 



We have used the same methods to calculate the smallest Carmichael number 
Sd with d prime factors for d up to 35. The results are given in Table|3| Heuristics 
suggest that a good approximation for \ogSd might be 21og2((i — l)log(d — 1). 
There is some support for this from Table 01 

In Table Eland Figure ^ we tabulate the function k{X), defined by Pomerance, 
Selfridge and Wagstaff ^ by 



C(X)^Xexp[ -k{X) 



^ log X log log log X 
log log X 

They proved that liminf A; > 1 and suggested that limsup/c might be 2, although 
they also observed that within the range of their tables k{X) is decreasing: Pomer- 
ance [HI, 13 gave a heuristic argument suggesting that lim/c = 1. The decrease in k 
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d 


N 


3 


561 


4 


41041 


5 


825265 


6 


321197185 


7 


5394826801 


8 


232250619601 


9 


9746347772161 


10 


1436697831295441 


11 


60977817398996785 


12 


7156857700403137441 


13 


1791562810662585767521 


14 


87674969936234821377601 


15 


6553130926752006031481761 


16 


1590231231043178376951698401 


17 


35237869211718889547310642241 


18 


32809426840359564991177172754241 


19 


2810864562635368426005268142616001 


20 


349407515342287435050603204719587201 


21 


125861887849639969847638681038680787361 


22 


12758106140074522771498516740500829830401 


23 


2333379336546216408131111533710540349903201 


24 


294571791067375389885907239089503408618560001 


25 


130912961974316767723865201454187955056178415601 


26 


13513093081489380840188651246675032067011140079201 


27 


7482895937713262392883306949172917048928068129206401 


28 


1320340354477450170682291329830138947225695029536281601 


29 


379382381447399527322618466130154668512652910714224209601 


30 


70416887142533176417390411931483993124120785701395296424001 


31 


2884167509593581480205474627684686008624483147814647841436801 


32 


4754868377601046732119933839981363081972014948522510826417784001 


33 


1334733877147062382486934807105197899496002201113849920496510541601 


34 


260849323075371835669784094383812120359260783810157225730623388382401 


35 


112505380450296606970338459629988782604252033209350010888227147338120001 



Table 3. The smallest Carmichael numbers with d prime factors 
for d up to 35. 



is reversed between 10^^ and 10^'*: see Figure^] We find no clear support from our 
computations for any conjecture on a hmiting value of k. 

In TableHwe also give the ratios C(10")/C(lO""i) investigated by Swift [Ul- 
Swift's ratio, again initially decreasing, also increases again before 10^'°^. 

In TableEland Figure[21we see that within the range of our computations C{X) 
is a slowly growing power of X: about X'^-^'^ for X around 10^^. 

In Table we give the number of Carmichael numbers in each class modulo m 
for m = 5, 7, 11 and 12. 

In Tables inland El we give the number of Carmichael numbers divisible by primes 
p up to 97. In Table |S| we count all Carmichael numbers divisible by p: in Table |5| 
we count only those for which p is the smallest prime factor. 

The largest prime factor of a Carmichael number up to 10^^ is p = 704988733, 
dividing n = 994018226608901845 = 5 • 13 • 1733 • 12517 • 704988733. We have 
p = n°-'^^^^: indeed n = p{2p - 1) and n - 1 = (p - l){2p + 1). 
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d 


log(6'd) 


21og2(rf-l)log(d-l) 


3 


3.293621187 


4 


2.324869052 


5 


1.772215418 


6 


1.755823184 


7 


1.503593258 


8 


1.385943139 


9 


1.296862582 


10 


1.273113126 


11 


1.210794211 


12 


1.187291825 


13 


1.183842403 


14 


1.142841324 


15 


1.115637251 


16 


1.112255293 


17 


1.068846695 


18 


1.086833722 


19 


1.067861876 


20 


1.055266652 


21 


1.056211783 


22 


1.041906608 


23 


1.034833298 


24 


1.024202006 


25 


1.026042173 


26 


1.014073506 


27 


1.017122489 


28 


1.010169037 


29 


1.007224804 


30 


1.000945160 


31 


0.984182030 


32 


0.993535534 


33 


0.990337138 


34 


0.984854273 


35 


0.984296312 



Table 4. The smallest Carmichael number with d prime factors 
Sd compared to 21og2((i — 1) \og{d — 1). 



The largest prime to occur as the smallest prime factor of a Carmichael number 
in this range is 891001, dividing 

804230935331967001 = 891001 • 895051 • 1008451. 

We note that this number is of the form (220fc+l)(221fc+l)(249fc+l) with k = 4050. 

We define the index of a Carmichael number N to be the integer i{N) — {N — 
1)/X{N) where A is the Carmichael function, the exponent of the multiplicative 
group. A result of Somer 'TU' implies that i{N) ^ cxd as iV runs through Carmichael 
numbers. In Table [Till we list the Carmichael numbers known to have index less than 
100: the list is complete for numbers up to 10^®. 

We define the Le/imer mdexof anumber A^ to be i'(A^) = {N—l)/(j){N). Lehmer's 
totient problem asks whether there is a number of integral Lehmer index greater 
than 1 (that is, whether there is a composite N such that (j){N) divides A^ — 1). Such 
a number would have to be a Carmichael number. See Guy |2j §B37 for further 
details. 
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Figure 1. k{X) versus X (expressed in powers of 10). 



n 


A: (10") 


C(10")/C(10"-^) 


3 


2.93319 




4 


2.19547 


7.000 


5 


2.07632 


2.286 


6 


1.97946 


2.688 


7 


1.93388 


2.441 


8 


1.90495 


2.429 


9 


1.87989 


2.533 


10 


1.86870 


2.396 


11 


1.86421 


2.330 


12 


1.86377 


2.286 


13 


1.86240 


2.339 


14 


1.86293 


2.319 


15 


1.86301 


2.353 


16 


1.86406 


2.335 


17 


1.86472 


2.373 


18 


1.86522 


2.394 



Table 5. 
n< 18. 



The function k{X) and growth of C{X) for X = 10", 



m 



The Carmichael numbers up to 10^^ with Lehmer index > 2 are Hsted in Table 
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n 


log C (10") /(n log 10) 


4 


.21127 


5 


.24082 


6 


.27224 


7 


.28874 


8 


.30082 


9 


.31225 


10 


.31895 


11 


.32336 


12 


.32633 


13 


.32962 


14 


.33217 


15 


.33480 


16 


.33700 


17 


.33926 


18 


.34148 



Table 6. C{X) as a power of X. 



m 


c 


25.10'-* 


10^^ 


10^^ 


lo^y 


10^4 


10^'^ 


10^" 


10^7 


5 





203 


312 


627 


1330 


2773 


5814 


12200 


25523 




1 


1652 


2785 


6575 


15755 


37467 


90167 


215713 


520990 




2 


82 


154 


327 


702 


1484 


3048 


6094 


12912 




3 


102 


172 


344 


725 


1463 


3059 


6285 


12845 




4 


124 


182 


368 


767 


1519 


3124 


6391 


13085 


7 





401 


634 


1334 


2774 


5891 


12691 


27550 


60602 




1 


1096 


1885 


4613 


11447 


28001 


69131 


168856 


414999 




2 


105 


186 


432 


967 


2109 


4599 


10011 


21961 




3 


152 


232 


496 


1055 


2178 


4707 


10039 


21918 




4 


129 


211 


450 


985 


2122 


4592 


9944 


21832 




5 


138 


222 


454 


1033 


2224 


4777 


10125 


21911 




6 


142 


235 


462 


1018 


2181 


4715 


10158 


22132 


11 





335 


547 


1324 


3006 


7032 


16563 


38576 


90731 




1 


640 


1131 


2770 


6786 


16548 


40891 


100071 


248840 




2 


139 


217 


473 


1068 


2361 


5338 


12054 


27268 




3 


142 


220 


457 


1045 


2348 


5319 


12186 


27483 




4 


104 


187 


442 


1026 


2317 


5261 


11917 


27276 




5 


152 


243 


466 


1066 


2370 


5316 


12194 


27670 




6 


116 


198 


440 


1061 


2400 


5384 


12155 


27362 




7 


122 


195 


458 


1023 


2223 


5165 


11853 


26960 




8 


129 


222 


475 


1107 


2450 


5449 


12012 


27639 




9 


131 


218 


465 


1042 


2285 


5179 


11835 


27158 




10 


153 


227 


471 


1049 


2372 


5347 


11830 


26968 


12 


1 


2071 


3462 


7969 


18761 


43760 


103428 


243382 


579215 




3 








1 


2 


2 


5 


5 


7 




5 


20 


32 


64 


124 


228 


448 


805 


1470 




7 


47 


75 


147 


289 


547 


1027 


1906 


3589 




9 


25 


36 


60 


103 


165 


294 


560 


1018 




11 














4 


10 


25 


56 



Table 7. The number of Carmichael numbers in each class modulo 
m ioT m — 5, 7, 11 and 12. 
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0.35 




Figure 2. C{X) as a power of X (expressed in powers of 10). 



p 


25.10^ 


10^^ 


10^^ 


10^^ 


10^^ 


10^'^ 


10^^ 


10^' 


10^^ 


3 


25 


36 


61 


105 


167 


299 


565 


1025 


1967 


5 


203 


312 


627 


1330 


2773 


5814 


12200 


25523 


53856 


7 


401 


634 


1334 


2774 


5891 


12691 


27550 


60602 


133380 


11 


335 


547 


1324 


3006 


7032 


16563 


38576 


90731 


214147 


13 


483 


807 


1784 


3998 


9045 


20758 


47785 


110320 


258539 


17 


293 


489 


1182 


2817 


6640 


16019 


38302 


91783 


221430 


19 


372 


608 


1355 


3345 


7797 


18638 


44389 


106273 


256020 


23 


113 


207 


507 


1282 


3135 


7716 


18867 


46612 


116538 


29 


194 


336 


832 


2094 


5158 


12721 


31110 


76647 


190486 


31 


335 


571 


1320 


3086 


7270 


17382 


41440 


99845 


242813 


37 


320 


535 


1270 


2926 


6826 


16220 


38647 


92744 


224582 


41 


227 


390 


1001 


2418 


5896 


14344 


34759 


84977 


209323 


43 


184 


296 


772 


1920 


4663 


11594 


28650 


70904 


176180 


47 


53 


80 


199 


492 


1223 


2873 


6810 


17002 


42440 


53 


92 


160 


351 


813 


2041 


5143 


12256 


30560 


76608 


59 


26 


41 


92 


262 


644 


1611 


3959 


9613 


23883 


61 


269 


453 


1075 


2542 


6047 


14429 


34503 


83268 


202835 


67 


110 


178 


407 


1063 


2540 


6306 


15295 


37521 


93475 


71 


104 


194 


521 


1320 


3351 


8546 


21485 


53857 


136174 


73 


198 


348 


849 


2145 


4925 


11929 


29072 


71137 


174025 


79 


64 


107 


247 


686 


1728 


4318 


10693 


26766 


66867 


83 


14 


24 


56 


137 


340 


838 


1929 


4640 


11465 


89 


68 


131 


320 


788 


1951 


4981 


12178 


30737 


77657 


97 


123 


193 


495 


1277 


3123 


7594 


18706 


45542 


112300 



Table 8. Primes occurring in Carmichael numbers. 
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p 


25.10'J 


10^^ 


10^^ 


10^-^ 


10^4 


10^'^ 


10^*^ 


10^' 


10^^ 


3 


25 


36 


61 


105 


167 


299 


565 


1025 


1967 


5 


202 


309 


624 


1325 


2765 


5797 


12175 


25481 


53776 


7 


364 


579 


1218 


2557 


5461 


11874 


25915 


57459 


127453 


11 


263 


428 


1071 


2509 


5979 


14397 


33893 


80745 


192905 


13 


237 


431 


1058 


2462 


5699 


13514 


32025 


76256 


183523 


17 


117 


206 


496 


1318 


3244 


8114 


20206 


50170 


125098 


19 


152 


244 


532 


1401 


3358 


8141 


20020 


49413 


122452 


23 


37 


78 


207 


535 


1360 


3317 


8195 


20803 


53549 


29 


55 


103 


284 


729 


1822 


4659 


11577 


29149 


73645 


31 


101 


168 


390 


876 


2116 


5153 


12575 


30667 


75541 


37 


60 


95 


219 


551 


1401 


3418 


8594 


21382 


52989 


41 


35 


68 


171 


414 


1092 


2736 


6788 


17275 


43376 


43 


35 


65 


168 


403 


943 


2308 


5520 


13636 


34030 


47 


14 


16 


36 


81 


195 


459 


1135 


2854 


7171 


53 


19 


30 


55 


147 


363 


973 


2327 


5842 


14731 


59 


2 


4 


11 


43 


100 


272 


618 


1542 


3825 


61 


34 


58 


148 


364 


851 


1978 


4722 


11278 


27910 


67 


8 


18 


50 


123 


317 


815 


1950 


4843 


12109 


71 


15 


25 


66 


161 


389 


979 


2480 


6178 


15602 


73 


14 


28 


68 


175 


406 


1015 


2508 


6277 


15612 


79 


4 


10 


17 


66 


175 


467 


1163 


2873 


7002 


83 


1 


1 


4 


8 


39 


79 


175 


457 


1143 


89 


10 


16 


23 


55 


148 


409 


1003 


2523 


6293 


97 


10 


20 


50 


106 


261 


606 


1413 


3445 


8448 



Table 9. Primes occurring as least prime factor in Carmichael numbers. 
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i 


N 


factors 


5 


6601 


7 • 23 • 41 


7 


561 


3- 11- 17 


18 


42018333841 


11 -47 -1049 -77477 


18 


55462177 


17 • 23 • 83 ■ 1709 


18 


8885251441 


11 -47 -1109 -15497 


21 


10585 


5 • 29 • 73 


22 


2465 


5 -17 -29 


23 


1105 


5 -IS- 17 


25 


11921001 


3 • 29 • 263 • 521 


31 


62745 


3 • 5 • 47 • 89 


37 


11972017 


43 • 433 • 643 


37 


67902031 


43 • 271 ■ 5827 


39 


334153 


19 • 43 • 409 


43 


52633 


7 • 73 • 103 


44 


15841 


7 • 31 • 73 


45 


8911 


7 -19 -67 


47 


2821 


7 -13 -31 


48 


1729 


7 -13 -19 


49 


1208361237478669 


53 • 653 • 26479 • 1318579 


50 


4199932801 


29 • 499 • 503 • 577 


52 


206955841 


17 -71 -277 -619 


53 


1271325841 


17 -31 •179- 13477 


54 


4169867689 


13 • 29 • 383 • 28879 


55 


271794601 


13 • 19 • 743 • 1481 


60 


6840001 


7 -17 -229 -251 


61 


1962804565 


5 • 103 • 149 • 25579 


67 


11985924995083901 


29 -101 -1427 -16349 -175403 


67 


410041 


41 - 73 • 137 


70 


162401 


17-41-233 


76 


4752717761 


11-17-107-173-1373 


81 


35575075809505 


5 - 197 - 223 - 353 • 458807 


82 


1496405933740345 


5 • 47 - 317 - 40253 - 499027 


83 


142159958924185 


5-37- 107-58379- 123017 


83 


158664761899885 


5 - 37 - 107 - 53987 - 148469 


83 


204370370140285 


5-37- 107-48677-212099 


83 


24831908105124205 


5 - 29 - 719 - 3023 - 78790717 


85 


171189355538562901 


23 - 71 - 983 - 1031 - 103437589 


90 


3778118040573702001 


11-47-1051-67967-102302009 


92 


520178982961 


11-29-131-607-20507 


94 


12782849065 


5 - 7 - 269 - 317 - 4283 


95 


8956911601 


11-17-127- 131-2879 


97 


5472940991761 


199 - 241 - 863 - 132233 


97 


721574219707441 


167 - 241 - 5039 - 3557977 


97 


83565865434172201 


103 - 1993 - 9551 - 42622169 


97 


9729822470631481 


127-409-110681-1692407 


99 


438253965870337 


43 - 139 - 49409 - 1484009 



Table 10. Carmichael numbers with small index 



10 
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I 


N 


factors 


2.00163 


3852971941960065 


3 


5 


23 


89 


113 -1409 -788129 


2.00388 


655510549443465 


3 


5 


23 


53 


389 • 2663 ■ 34607 


2.00438 


978502408508476545 


3 


5 


23 


89 


173 ■ 197 • 13553 • 68993 


2.00837 


13462627333098945 


3 


5 


23 


53 


197 ■ 8009 • 466649 


2.01001 


26708253318968145 


3 


5 


17 


lis 


• 57839 • 16025297 


2.01548 


121719617715279585 


3 


5 


17 


89 


5897 -6833 -133103 


2.02272 


371717470696746945 


3 


5 


29 


47 


113 -197 -353 -449 -5153 


2.03067 


269040992399565 


3 


5 


23 


29 


4637 ■ 5799149 


2.03207 


158353658932305 


3 


5 


17 


89 


149 • 563 • 83177 


2.03614 


1817671359979245 


3 


5 


23 


29 


359 -11027 -45893 


2.05825 


184882434303977985 


3 


5 


17 


47 


113 - 449 • 304041473 


2.07138 


16057190782234785 


3 


5 


17 


29 


269 - 6089 - 1325663 


2.07294 


75131642415974145 


3 


5 


23 


29 


53 -617 -9857 -23297 


2.08045 


202703674932642885 


3 


5 


23 


29 


53 • 197 • 3389 - 572573 


2.08129 


881715504450705 


3 


5 


17 


47 


89 - 113 - 503 - 14543 


2.08283 


462937246809774465 


3 


5 


17 


29 


107-65729-8901089 


2.08287 


781506906921758865 


3 


5 


17 


23 


32117-57077-72689 


2.08309 


224811969542371905 


3 


5 


17 


29 


113-2297-6977-16787 


2.08381 


31454143858820145 


3 


5 


17 


23 


2129 - 39293 - 64109 


2.08894 


6128613921672705 


3 


5 


17 


23 


353 - 7673 - 385793 


2.09063 


220629545178715905 


3 


5 


17 


23 


353 - 1409 - 7673 - 9857 


2.09805 


947087538769733505 


3 


5 


17 


23 


173 - 929 - 2237 - 449153 


2.09841 


12301576752408945 


3 


5 


23 


29 


53-113- 197-1042133 


2.11432 


1886616373665 


3 


5 


17 


23 


83 - 353 - 10979 


2.11493 


180950795673242145 


3 


5 


17 


23 


83 - 353 - 2663 - 395429 


2.11493 


3193231538989185 


3 


5 


17 


23 


113-167-2927-9857 


2.12088 


916541202603455265 


3 


5 


17 


23 


89 - 347 - 353 - 947 - 15137 


2.13936 


11947816523586945 


3 


5 


17 


23 


89-113-233-617-1409 


2.14083 


101817952350880305 


3 


5 


17 


23 


89-113-149-3257-3557 


2.23494 


171800042106877185 


3 


5 


17 


23 


29 - 53 • 89 - 197 - 1086989 



Table 11. Carmichacl numbers with Lehmer index greater than 2 
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